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We consider soft-core bosons with onsite interaction loaded in the honeycomb lattice with different 
site energies for the two sublattices. Using both a mean-field approach and quantum Monte-Carlo 
simulations, we show that the topology of the honeycomb lattice results in a non-vanishing Berry 
curvature for the band structure of the single-particle excitations of the system. This Berry curvature 
induces an anomalous Hall effect. It is seen by studying the time evolution of a wavepacket, namely 
a superfluid ground state in a harmonic trap, subjected either to a constant force (Bloch oscillations) 
or to a sudden shift of the trap center. 


I. INTRODUCTION 

Topology and gauge fields, exemplified by Chern num¬ 
bers, Berry phases, Zak phase and Berry curvatures m, 
are key concepts at the heart of many condensed-matter 
phenomena mg. For solid-state systems, Bloch’s theo¬ 
rem introduces wave numbers k belonging to a parameter 
space with the topology of a torus, the Brillouin zone, 
and a set of periodic wave functions depending para¬ 
metrically on k which offer natural settings for bundles 
and connections 0. Recently, ultracold atomic systems 
have progressively and successfully confirmed their abil¬ 
ity to mimic or emulate some paradigmatic phenomena 
of condensed-matter systems, in particular topological ef¬ 
fects. Indeed one can now load independent or interact¬ 
ing bosons and/or fermions into two-dimensional, very 
well controlled, optical lattices [8HI2] and one can gener¬ 
ate carefully designed synthetic magnetic fields [EM]. 
These advances pave the way to accurate cold atoms 
experiments targeting topology-related effects such as 
defects m, color superfluidity [18], momentum-space 
Berry curvatures [19] or quantum Hall states with strong 
effective magnetic fields [2(11 [25] . 

Because of its remarkable low-energy electronic exci¬ 
tations, graphene has been the source of many key dis¬ 
coveries [26] (27 which sparked a vivid research flow now 
reaching new territories as exemplified by ultracold atoms 
loaded in optical lattices [24l [281(361 . In particular, the 
tight-binding model on the honeycomb lattice or equiv¬ 
alently on the brick wall lattice is well known to exhibit 
a Berry curvature in its band structure oisaisH], such 
that the anomalous Hall effect can be observed 0[39]. In 
this paper, we discuss the impact on the Berry curvature 
of the interactions between bosons loaded in the honey¬ 
comb lattice with site energy imbalance, i.e. when the 
site energies of the two sub-lattices are different. In that 
case, the phase diagram depicts Mott-insulating phases 
at half-integer filling, a feature that is shared with the 


energy imbalanced square lattice. On the other hand, as 
we show below, the properties of the Bogoliubov excita¬ 
tions above the ground state are different for these two 
lattices: only the excitations for the honeycomb lattice 
depict a non-vanishing Berry curvature. This is empha¬ 
sized by looking at a dynamical situation, such as the 
Bloch oscillations, where the anomalous Hall effect is ob¬ 
served. 

The paper is organized as follows. In section [ilj we 
introduce the model, summarize the basic properties of 
the energy imbalanced honeycomb lattice, in particular 
the Berry curvature and we present the two approaches 
used in the paper: the quantum Monte-Carlo method 
(QMC) [4014431 and the Gutzwiller ansatz (GA) [441446] . 
In section m we present our numerical results: the 
ground state phase diagram and the excitations, empha¬ 
sizing that they depict non-vanishing Berry curvature. 
In section [TV| we explain how to observe the anomalous 
Hall effect resulting from the Berry curvature in the exci¬ 
tations. A summary of results and conclusions are given 
in section El 

II. MODEL 

A. Single-particle Hamiltonian and Berry 
Curvature 

We start with a general bipartite lattice structure with 
four nearest neighbors. Each primitive unit cell con¬ 
tains two lattice sites A and B , as depicted in Fig. [l] 
The whole lattice structure can be generated by repeated 
translations of a unit cell along the Bravais primitive vec¬ 
tors 

o,\ — €> x and Q >2 — &x H - &y (l) 

where we have assumed the spacing between two sites in 
the primitive unit cell is unity. The primitive reciprocal 
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lattice vectors are 


and eigenfunctions [38] 


i>i = 7r(e x - e y ) and b 2 = 7r(e x + e y ), ( 2 ) 

and fulfills the relation cq • bj = 2iT§ij. 

The displacements that move a A site to a neighboring 
B site are parameterized by bond vectors 
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FIG. 1. (Color online) Structure of a bipartite lattice with 
four nearest neighbors. Sites belonging to sublattice A are 
represented by dark (blue) points in the figure, whereas B 
sites are represented by light grey (green) points. The hopping 
amplitudes are all equal in (a) while in (b) the hopping term 
£4 along links C4 is different from the other three near neighbor 
hopping terms. 


Atoms trapped in the lattice can hop from one site 
to the neighboring sites with amplitudes t v [y = 1..4) 
corresponding to the four bond vectors. In the absence 
of interactions, the Hamiltonian in real space reads 


H {) = - ^ (tij b\bj + t*j bjbij - n i 

( i,j) iEA j£B 

( 4 ) 

where E {£ 1 , £ 2 , £3, £4}, bi is the operator annihilat¬ 
ing one particle on lattice site i, rii = b\b i is the atom 
number operator, jiA and fis are chemical potentials for 
sublattice A and B respectively. When all hopping am¬ 
plitudes t v are the same, Eq. @ describes a standard 
square lattice. If £4 = 0, the lattice becomes the brick 
wall lattice which is topologically equivalent to the hon¬ 
eycomb lattice. Using Fourier transform, we obtain the 
Hamiltonian in momentum space 
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where k is the Bloch wave vector, A = — (/jla ~ Hb)/ 2, 
M = (ma + Ms)/2, N = Ek( & b & kA + 4s & ks) is the 
total number of particles, and Tk = — Ylt=i ^ e_2kC ” • 
The Hamiltonian ([5| has the eigenvalues 


(6) 


where (/) = — argTk and /? = arcsin( |Tk |/\/^ 2 + |Tk | 2 )• 
For A = 0, the gap between the two energy bands closes 
at |Tk| = 0, and exhibit the so-called Dirac cone. To open 
a gap, one needs to imbalance the chemical potential be¬ 
tween the A and B sites. In the following, we always 
consider the gapped case, i.e. A^O. 

The Berry curvature for each band of the Hamilto¬ 
nian can be calculated as [38] 


n±(k) = ±\ 
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Expanding the Bloch wave vector k over the reciprocal 
primitive vectors, k = aqiq + a 2 b 2 , we find the Berry 
curvature can be written as follows 


Ofc(k)«±I--- m 

2 (A 2 + |r k | 2 ) 3/2 


x [tit3 — t^ti] sin(27ro:2) — [ii ^2 — Ms ^ 4 ] sin (27 tq;i) 

( 9 ) 

where we have assumed real values for the hopping am¬ 
plitudes. For the square lattice (t u = £), the Berry cur¬ 
vature is always zero, reflecting the presence of both the 
inversion and the time-reversal symmetry in the system. 
However, if one of the hopping amplitudes is different 
from the other three, for example £4 7^ £ 1 , 2 , 3 , the sym¬ 
metry under the exchange x -a —x is broken. As a con¬ 
sequence, the system can exhibit a non-vanishing Berry 
curvature. In this case, as the symmetry under the ex¬ 
change y -A —y remains, /3( k) and cj)(k) are even func¬ 
tions of k y . According to Eq. pi), one can easily see that 
the Berry curvature should be an odd function of k y , i.e. 
(kx^ky^) — ky) • 

In Fig. [2| we give an example of the Berry curvature 
in the first Brillouin zone (BZ), i.e., for |c^| < 1/2, calcu¬ 
lated from § for the imbalanced honeycomb lattice, i.e. 
£ 1 , 2,3 =£,£4 = 0 and A = 5£. Due to the time-reversal 
symmetry, the Berry curvature has the usual property 
fi±(—k) = —fi±(k) such that, for each band, the Chern 
number which is the surface integral of Berry curvature, 
7± = f BZ dSft±(k), is vanishing. This emphasizes that 
the bands are topologically trivial. Furthermore, due to 
the odd symmetry of the Berry curvature under the ex¬ 
change k y -A —ky , it is also an even function of k x , 
£l±(k x ,k y ) = Q±(—k x ,k y ), as shown in Fig. |2j in the 
(ai,a 2 ) plane, the Berry curvature is symmetric with 
respect to the anti-diagonal and antisymmetric with re¬ 
spect to the diagonal. 


e ±(k) = ±\/ A 2 + |r k | 2 — ji 
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FIG. 2. (Color online) Berry curvature in the first Brillouin 
zone for the honeycomb lattice calculated from Eq. § . The 
parameter are chosen as follows: £1,2,3 = £, £4 = 0, A = 5£. 


B. Numerical methods 

When the interaction on a site is taken into account, 
the system can be described by the Bose-Hubbard Hamil¬ 
tonian 


= + (10) 

i 

where Hq is given by Q , and U is the on-site interaction 
strength. This is the Hamiltonian we shall investigate 
in this study using two complementary methods that are 
described below. 


1. Quantum Monte-Carlo 

The ground state of Hamiltonian ( fl0| ) can be stud¬ 
ied using the stochastic series expansion (SSE) quan¬ 
tum Monte Carlo method with operator-loop updates 
[4QH42] . The SSE is a finite-temperature QMC algo¬ 
rithm based on importance sampling of the diagonal ma¬ 
trix elements of the density matrix e~^ n . There are 
no approximations beyond statistical errors. Using the 
operator-loop cluster update, the autocorrelation time 
for the system sizes we consider here is at most a few 
Monte Carlo sweeps for the entire range of parameter 
space explored. The simulations are carried out on finite 
lattices with L 2 sites for L up to 32 at temperatures suffi¬ 
ciently low in order to resolve ground-state properties of 
this finite system [43] . Estimates of physical observables 
in the thermodynamic limit are obtained from simulta¬ 
neous finite-size and finite-temperature extrapolation to 
the L —)> 00 , /3 —)> 00 limit. We consider here the soft-core 
boson case which allows multiple occupation to occur (up 
to n( max ) = 4 bosons per lattice site are allowed in the 
present study). Periodic boundary conditions along the 
x and y axis have been used. 


2. Mean-field 

A well known mean-field method to solve the Bose- 
Hubbard model is the Gutzwiller ansatz [471450] . where 
the ground state wavefunction is assumed to be a tensor 
product of on-site wavefunctions: 


N 

W = 0 I ipj) wh ere | ipj) = ^ fnj\n,j) (11) 

j n= 0 

where | n,j) represents the Fock state of n atoms occu¬ 
pying the site j, N is a cutoff in the maximum number 
of atoms per site, and f n j is the probability amplitude 
of having the site j occupied by n atoms. 

Minimizing the mean-field energy (\k|i4|\h) over the 
amplitudes f n j allows us to determine the mean-field 
ground state properties as functions of the different pa¬ 
rameters ([/, £, /ia, Hb)- For instance, the superffuid 
phase corresponds to a non-vanishing value of the order 
parameter (\h|^|\h), whereas the Mott phase corresponds 
to a vanishing order parameter and the | ipi) are pure Fock 
states. 

In addition, the preceding ansatz can be extended 
to the time domain giving us access not only to the Bo- 
goliubov excitations above the ground state, but also to 
the full mean-field evolution of the interacting wavefunc¬ 
tion. More precisely, the time evolution is obtained by 
solving the following set of equations: 

d t df* d ■ ( } 

The Bogoliubov excitations are obtained by expanding 
the amplitudes / n j(£) around their ground state values 

fnf namely: 


fn,j ~ fnj + 9n,j(t ) 




(13) 


where ujj is the frequency of the ground state evolution. 
Assuming \g n ,j{t)\ \fn,j\ an d keeping only the linear 
terms in the dynamical Eqs. (12), one obtains the usual 
Bogoliubov equations: 
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(14) 


where g(t) is a shorthand notation for the vector 
(• • • , go,j 5 gij, ••• , gN,j 7 • • • ) T and £ Fas the usual Bo¬ 
goliubov structure: 


£ = 




(15) 


A and B are complex matrices satisfying = A and 

B T = B. 

Due to the U( 1) invariance, the values of can be 
taken real. Furthermore, as we have observed (see be¬ 
low), the mean-field ground state does not break the 
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translation invariance, such that the values are the 
same for all equivalent sites, namely /^°] = A n for all 


A —sites and /^°j = B n for all B— sites. Because of trans¬ 
lation invariance, the Bogoliubov Eqs. (14) can be di¬ 
agonalized in momentum space, taking into account the 
bipartite nature of the lattice: 


9n,j(t) = ^V (k '^ Wt)uk >A' 

k 

9n,i(t) = E ei(k ' B '' Wt)uk - S ^ 


+ e -i^A 3 -.t )<An{w) 


for A sites and B sites respectively, leading to the follow¬ 
ing eigensystem 


< i8 > 

where u k = (^k,A,o, * • • , ^k,A,;v, ^k,B,o, * • • , ^k ,b,n) and 
^k = (Dc,A, 0 , • • • 5 Vk,A,N•> ^k,B,0> " * , ^k ,B,n)- £k has the 
same structure as £, which leads to the usual property: if 
(iq v) is an eigenstate for the energy a;, then (v*,u*) is an 
eigenstate for the energy —uj*. Therefore one can focus 
on the eigenstates with positive skew-norm only, i.e. such 
that u^u — v^v = 1. The ground state is stable if all these 
eigenstates have real and positive eigenenergies co>k- 


III. RESULTS 
A. Phase diagram 

Figure [3] shows the ground state phase diagram for 
a honeycomb lattice (£4 = 0 ) obtained from ansatz 
Eq.jn). Due to different chemical potentials between 
the A and B sites, the insulating lobes occur not only 
at (n) = 0 , 1, 2 , ... but also at (n) = 1 / 2 , 3/2, — We 
compare the mean-field result with QMC simulation in 
Fig- i where we plot the density per site as a func¬ 
tion of fi = (ha + I^b)/ 2 , fixing the hopping amplitude 
at t/U = 0.03, which corresponds to the white lines in 
Fig-! Very good agreement is found between the mean- 
field approach and QMC simulations. 


B. Excitations 

The phase diagram in Fig. [3] itself is not a signature 
of the Berry curvature induced by the topology of the 
honeycomb lattice. Indeed, the Hubbard model on the 
square lattice (i.e. with a non-vanishing hopping param¬ 
eter £4) with the same imbalance /a a — Hb exhibits a 
similar phase diagram, including the half-integer filling 
Mott phases. The reason is that the ground state corre¬ 
sponds, roughly speaking, to k = 0, precisely where the 
Berry curvature D(k) vanishes. Therefore, the impact 
of the Berry curvature has to be found in the excitation 
properties. 


(a) 
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FIG. 3. (Color online) Mean-field phase diagram of the Bose- 
Hubbard Hamiltonian Eq. (10) on the honeycomb lattice, 
(a) Density per site; (b) the superfluid order parameter. The 
parameters used in the calculation are £1,2,3 = £, £4 = 0, 
A/U = 0.15, N = 10. The white lines correspond to t/U = 
0.03. 


As explained above, within the mean-field approach, 
the Bogoliubov excitations are obtained from the diag- 
onalization of £k, leading to a band structure for the 
Bogoliubov spectrum. For instance, in Fig. [5j we show 
the lowest excitation band for the honeycomb lattice for 
two different phases, namely the Mott phase (/jl/U = 0 ) 
and the superfluid phase (/jl/U = 0.15). The rest of the 
parameters are the same as in Fig. [4j As expected, the 
excitations are gapped in the Mott phase, whereas, in 
the superfluid phase, they are gapless and exhibit a lin¬ 
ear behavior at small momenta. 

In general the Bogoliubov bands are isolated allowing 
the numerical computation m of the Berry curvature 
associated with each band, as shown in Fig. [ 6 j for the two 
excitation bands of Fig. [5j One clearly sees that in both 
situations, the Berry curvature is non-vanishing, being 
maximum (in absolute value) along the anti-diagonal, i.e. 
the k y axis, like in the non-interacting case. In addition, 
the time reversal symmetry and the y —>> —y symmetry, 
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FIG. 4. (Color online) Comparison of QMC and mean-field 
results for the density per site as a function of total chemical 
potential /i — (fi A + I^b)/ 2. Parameters are £1,2,3 = £, £4 = 
0, A/U = 0.15, as in Fig. 3] and the hopping amplitude is 
fixed to t/U = 0.03. The blue solid curve is obtained from 
the mean-field approach (corresponding to the white line in 
Fig.§, the red dots are obtained from QMC simulations. The 
lattice size for QMC simulations is 16 x 16, with a maximum 
atom number per lattice site n^ max ^ = 4. 


like in the non-interacting case, imply that the Berry 
curvature is odd under both transformations k —>> — k and 
k y —>> —k y and that it is even under the transformation 
k x —> —k x . In both Fig. [6] and Fig. [5| these properties 
are clearly seen: (i) antisymmetric under inversion and 
with respect to the diagonal; (ii) symmetric with respect 
to the anti-diagonal. 

Similarly, one could define a pseudo Berry curvature 
for the excitations obtained from QMC computations: 
because of the bipartite nature of the honeycomb lattice, 
the equal time Green function in k-space is a 2 x 2 matrix: 


G(k) = 


G AA (k) G ab ( k)\ 
GW(k) Gbb{ k) J ’ 


(19) 


where Gll'Qz) = (b^ L bw)- Since G(k) is hermitian, 
one can write G(k) = go(k)t + g(k) • cr. where a = 
( a x , cr y ,cr z ) are the Pauli matrices. Thereby, one can de¬ 
fine the spherical angles /3 g ( k) and cf) g ( k), corresponding 
to the direction of g(k), or, equivalently, obtained from 
the diagonalization of G(k), similarly to diagonalizing 
Hq, see Eq.([5| and Eq.([7|). From these two angles, one 
can define a pseudo Berry curvature: 


«p(k) = o 


1 ( d cos /3 g d(j)g d cos /3 g d(j) Q 


dk r dk v 


dk y dk x 


( 20 ) 


This curvature, being computed from the equal-time 
Green function, is not directly related to a genuine Berry 
phase that could be measured from the excitations of the 
system. For that, one would need to define a curvature 
from the real frequency Green function G(k, ui). That en¬ 
tails computing the imaginary time Green function and 
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FIG. 5. (Color online) Bogoliubov excitation spectrum u;(k). 
Parameters are the same as in Fig. [4] The top panel cor¬ 
responds to the Mott Phase g/U = 0 and the bottom one 
to the superfluid phase g/U = 0.15. As expected, the excita¬ 
tions are gapped in the Mott phase, whereas, in the superfluid 
phase, they are gapless and exhibit a linear behavior at small 
momenta. 




subsequent analytic continuation to the real axis, which 
is well known to be a delicate task. Therefore, in order 
to emphasize the impact of the lattice topology on the 
excitations, we shall focus on the pseudo Berry curvature 
defined in Eq. (20). The results are shown in figure [ 7 ] cor¬ 
responding to a ground state in the n = 1/2 Mott phase. 
For comparison, we also show the pseudo Berry curvature 
obtained in the same ground state but for the square lat¬ 
tice. Only the curvature for the honeycomb lattice is non¬ 
vanishing and is largest (in absolute value) around the 
same position as for the non-interacting case, i.e. around 
the location of the conical intersections. In addition, this 
pseudo Berry curvature has the same symmetry proper¬ 
ties: Q^(k) = — Q g (— k), (time reversal symmetry), odd 
under k y —>• — k y and even under k x —k x . 


In the superfluid phase, the 2x2 matrix is dominated 
by its diagonal elements and the resulting pseudo Berry 
curvature is smaller than the QMC accuracy. 
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FIG. 6. (Color online) Berry curvature of the Bogoliubov 
excitations displayed in Fig. [5] One can clearly observe 
Q(k) = — H(—k) (time reversal symmetry), the antisymmetry 
with respect to the diagonal (k y —>• — k y ) and the symmetry 
with respect to the anti-diagonal (k x —>• —k x ). 


FIG. 7. (Color online) Pseudo Berry curvature, Eq.(20), ob¬ 
tained with the QMC equal-time Green functions for the hon¬ 
eycomb and square lattices. The upper panel corresponds 
to the n — 1/2 Mott phase for the honeycomb lattice, i.e. 
£1,2,3 = £, £4 =s 0, A/U — 0.15, ii/U — 0 and t/U — 0.03. 
The lower panel corresponds to the n — 1/2 Mott phase for 
the square lattice, i.e. £ 1 , 2 , 3,4 = £. For the square lattice, one 
obtains only a random pattern with an amplitude two orders 
of magnitude lower than the honeycomb lattice results. This 
emphasizes that the interacting system must really depict a 
non-vanishing Berry curvature in the excitation spectrum. In 
addition, one can see that the pseudo curvature has the same 
properties as the non-interacting system: it attains the largest 
(absolute) values around the location of the conical intersec¬ 
tions and fulfills £l g (k) = — Q g (— k), due to the time reversal 


symmetry. It is also odd under k y 

k x y k x . 


—k v and even under 
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IV. ANOMALOUS HALL EFFECT 

Although the preceding analysis has shown that the 
impact of the Berry curvature can be found in the prop¬ 
erties of the excitations, one can still raise the question 
whether it can be observed. It is well known that in this 
situation, the system is expected to exhibit an anoma¬ 
lous Hall effect (AHE) - a Hall effect without applying 
an external magnetic field - in the non-interacting case. 
In this section, we shall demonstrate that one can still 
observe the AHE in the presence of the interaction. How¬ 
ever, the AHE can be observed only if the initial state 
is a wavepacket well-localized in real space. Indeed, if 
one starts with a pure Bloch state, the effect of a con¬ 
stant force only leads to Bloch oscillations, i.e. a periodic 
variation of the quasi-momentum along the direction of 
the force, and the Berry curvature simply results in a 



modification of the phase accumulated along the path. 

A simple way to prepare a wavepacket is to start from 
the ground state in a harmonic trap, as depicted in Fig.[8j 
the atoms are well localized in space, still, the wavepacket 
is also well localized in the Brillouin zone around k = 0. 
The harmonic trap is such that the chemical potential at 
the center corresponds to the superfluid phase described 
in the previous section, i.e. y/U = 0.15 and t/U = 0.03. 
The site energy of the A-sublattice is lower than that of 
the H-sublattice, corresponding to A /U = —0.15, result¬ 
ing, in the present case, in a Mott-like state for the A-site 
with a flat density at unit filling, see Figj8j top-left. On 
the contrary, the filling of the F-sublattice is low, such 
that the system is in a superfluid phase, with a smooth, 
gaussian like, density, see FigjHJ top-right. In both cases, 
the density in the k-space is peaked at k = 0. 



FIG. 8. (Color online) Initial wavepacket in a harmonic trap. The atomic density in the A-sublattice (resp. B-sublattice) is 
displayed in the top left (resp. top right) plot. The harmonic trap is such that the chemical potential at the center corresponds 
to the superfluid phase described in the previous section. Since the chemical potential yA is larger than / ib , i.e. A/U = —0.15, 
the A-sites are in a Mott-like state with unit filling. On the contrary, the filling of the B-sublattice is low, such that the system 
is in a superfluid phase. Densities in the k-space for both the sublattices are displayed in the bottom row and are found to be 
peaked around k = 0. 


A. Bloch oscillations 


To simulate Bloch oscillations, starting from the 
wavepacket, one removes abruptly the harmonic trap and 
adiabatically increases the linear potential along the y 


axis at time r = 0 mimicking the effect of an electric field 
E along the y axis. The evolution of the wavepacket is 
obtained by solving the time-dependent equations for the 
mean-field amplitudes Eq.( 12). In k-space, this leads to a 
smooth evolution of the wavepacket (with some broaden- 







ing) along k y , see Fig. 9] This is emphasized in Fig. [lOj 
where we have plotted the average value of k y (within 
the H-sublattice) as as function of time. After the adi¬ 
abatic transition, the linear increase of k y with time is 
a clear signature of the Bloch oscillation. The AHE is 
demonstrated by the bottom plot in Fig. [lOj where we 
show the displacement from the center of the trap of the 
wavepacket along the x-axis, i.e. perpendicular to the ap¬ 
plied force. One can clearly see the impact of the Berry 
curvature, in particular when comparing with a similar 
evolution but for the square lattice, for which there is 
no displacement. In addition, changing the sign of the 
linear potential i.e. the direction of the applied force, re¬ 
sults in changing the direction of the Bloch oscillations 
k(r) —)> —k(r), i.e. changing the direction of the group 
velocity of the wavepacket. Together with the fact that 
the Berry curvature is odd with respect to the change 
k —>> — k, this results in a Hall displacement in the same 
direction along the x-axis, as depicted in Fig. [lOj the 
two curves coincide with each other. Note that we have 
checked that this effect is independent of the exact lo¬ 
cation of the center of the trap. Finally, we have also 
verified that if the force is applied along the x-axis, no 
AHE is observed since the Berry curvature vanishes for 
k y = 0. 





FIG. 10. (Color online) Bloch oscillations and AHE for a 
wavepacket. The top plot shows the average value of k y 
(within the B-sublattice) as a function of time. After the 
adiabatic transition, the linear increase of k y with time is a 
clear signature of the Bloch oscillation. The two curves, i.e, 
the continuous (black) line and the dashed (red) line with 
the circles, correspond to opposite directions of the effective 
electric field. The AHE is demonstrated by the bottom plot, 
where one plots the displacement from the center of the trap 
of the wavepacket along the x-axis, i.e. perpendicular to the 
applied force: the continuous (black) line with the squares and 
the dashed (red) line with the circles correspond to opposite 
directions of the effective electric field. For comparison, the 
corresponding displacement on a square lattice is also shown 
by the continuous (green) line with the diamonds. 



FIG. 9. (Color online) Evolution of the wavepacket in momen¬ 
tum space. The upper panel displays the wavepacket after a 
time r = 25/1/, the lower one after r = 50 /U. The center 
of the wavepacket moves at a constant velocity along the k y 
axis, i.e., an + 02 — 0. The broadening of the wavepacket is 
weak. 


B. Shift of the trap center 

The amplitude of the AHE discussed above is rather 
small. One way to get a larger (measurable) effect con¬ 
sists of shifting abruptly the center of the trap at r = 0 
along the y- axis, such that the wavepacket experiences a 
net force in that direction. Figure [ll] shows the evolution 
of the expectation value of the center of the wavepacket, 
i.e., (r), after a shift of the trap center by 20 lattice spac- 
ings in the y direction. As one can see, not only does the 
wavepacket evolve along the shifted center, but, at the 
same time, it moves along the transverse direction, i.e. 
the x-axis. Similarly to the Bloch oscillations, one has 
the following results: (i) nothing happens on the square 
lattice; (ii) nothing happens when the shift is along the x- 
axis where the Berry curvature always vanishes; (iii) the 
displacement along x is in the same direction whether the 
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trap center is shifted towards the +y or the —y direction. 



FIG. 11. (Color online) Evolution of the average center of 
the wavepacket, i.e., (r), after a shift of the trap center by 
20 lattice spacings in the y direction (continuous (black) line 
with the squares. The wavepacket clearly exhibits motion 
along the transverse direction, i.e. the rr-axis. Similarly to 
Bloch oscillations, one has the following results: (i) nothing 
happens on the square lattice, see the (green) line with the 
diamonds; (ii) nothing happens when the shift is along the x- 
axis where the Berry curvature always vanishes, see the (blue) 
line with the triangles; (iii) the displacement along x is in the 
same direction whether the trap center is shifted towards the 
-\-y or the —y direction, see the dashed (red) line with the 
circles. 


V. CONCLUSION 


In summary, we have shown that for bosons in the 
honeycomb lattice with energy imbalance, the Berry cur¬ 
vature originally seen in the band structure, i.e. for the 
non-interacting system, is also present in the excitations 
above the ground state of the interacting system, in both 
the Mott-insulating phase and the superfluid phase. In 
addition, we have shown that one consequence of the 
Berry phase, the anomalous Hall effect, could be observed 
in dynamical experiments, such as Bloch oscillations. 
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